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TECBMICAL NOTE 3701 

COMPARISON OP SEVERAL METHODS FOR OBTAINING THE TIME 
RESPONSE OP LINEAR SYSTEMS TO EITHER A UNIT IMPULSE OR 
ARBITRARY INPUT PROM EREQUENDY-RESFONSE DATA 
By James J. Donegan and Carl R. Huss 


SUMMARY 


Several methods of obtaining the time response of linear systems to 
either a unit Impulse or an arbitrary Input from frequency-response data 
are described and con:q)ared. 

Comparisons Indicate that an the methods give good accuracy vhen 
applied to a second-order system; the main difference Is the required 
conq)uting time. Several of the methods when applied to hi^er order 
systems require excessive computing time in order to obtain the same 
degree of accuracy. The methods generally classified as inverse Laplace 
transform methods were found to be most effective In determining the 
response to a \anit lirqjulse from frequency-response data of higher order 
systems . 

Some discussion and exaaq>les are given of the use of the methods 
as fli^t-data-analysis techniques In predicting loads and motions of a 
flexible aircraft on the basis of simple calculations when the aircraft 
frequency response is Imown. 


INTRODUCTION 


The frequency-response type of analysis used on linear systems has 
found extensive application in the field of aircraft stability and in 
the determination of overall dynamic characteristics of an aircraft. In 
faet the current trend to perform analysis on flight test data in frequency- 
response form appears to be gaining favor especially in the case of flex- 
ible aircraft. Considerable emphasis has, therefore, been given to methods 
of determining the frequency response of a system from transient responses. 
A concise re'sume and comparison of methods for obtaining the frequency 
response from transient responses are presented in reference 1. 

This paper briefly considers the next step In the process - the 
converting of the Information contained in a frequency response to the 
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time plane in the form of the response to a unit impulse. For certain 
purposes this conversion gives data in more useful form. The re^onse 
of a linear system to a unit tn 5 >ulse may be used in conjunction with 
Duhamel's (superposition) Integral to determine (l) the aircraft transient 
response to any type of input or (2) the Input required to cause any 
required, aircraft transient response. 

It appears that methods of converting frequency-response data to 
transient data present a f li^t-data-analy s is technique which permits 
the prediction of aircraft motions and loads for a flexible aircraft 
without knowledge of the equations of motion relating the input and 
oxxtput. These methods also bypass the need for computing transfer- 
function coefficients or stability derivatives in predicting these loads 
and motions. Such predictions are Inportant in anticipating the motions 
and loads for more hazardous aircraft maneuvers. These methods are also 
Tiseful In predicting time responses of complicated linear systems whose 
frequency response is known. 

The pmpose of this report is to collect and briefly compare a few 
of the methods now available for performing this operation. The methods 
are coopared on the basis of accuracy, computing time required, and 
applicability of the method to hi^er order systems. Some extensions of 
these methods are also given. With such information available, engineers 
may 'then select the method which best fits their needs . 


SYMBOLS 


F(t) 

function of time 


H(icu) 

frequency response of a system 


H(s) 

transfer function (in terms of 
relating input and output 

Laplace variable 

h(t) 

time response of a system to a 

unit impulse 

II 

•H 



J 

index of summation 



transfer coefficients defined by equation (8) 


K,L limits on summation of P-transform 


M 


Mach nxmiber 
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m 

index of simmtion of P-ta^ansf orm. 

N 

•botal number of da*ba points used in solution of 
equation (9) 

n 

limit of summation 

P 

polynominal transform operator 

Re |H(ia))] 

real part of frequency response defined "by equation (21) 

r 

ordinate of pxolse used to fit RejnCicoj] 

s 

Laplace variable, c + ico 

t 

time, sec 

X 

output 

x(s) 

Laplace transform of x(t) 

x(t) 

response of a system to an arbitrary input 

1 ^1 
1 6| 

anqolitxide ratio of frequency response 

A 

Indicates increment 

5 

input 

8(t) 

time history of liq)ut to a linear system 

d 

pitching velocity, radians/sec 

T 

value of time, sec 

i 

phase angle, deg 

^x8 

phase angle between output x and input 6 of frequency 
response; negative |hase angles indicate lag 


clrcxilar frequency, radians/sec 

“d 

damped natural frequency, radians/sec 


fundamental frequency, radians/sec 


undamped natural frequency, radians/sec 
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Subscripts : 
e elevator 

1 index of summation indicating rov 

J index of summation indicating column 

^ A tilde indicates polynomlnal transform of function; for example, 
F(x) denotes polynomlnal transform of F(t). 

DESCRIPTION AND D3SCUSSI0N OF MEEBDDS 


In this section the methods far determining the response to a unit 
impulse aaad the methods for obtaining the response to arbitrao^y inputs 
are discussed. 

The methods for obtaining the time response to a unit impulse from 
frequency-response data fall into several basic categories -which may be 
gener ally classified as follows; 

( 1 ) Inverse Laplace transform methods 

( 2 ) Fourier method 

(3) Other methods 

These methods are briefly outlined to indicate the techniques involved.' 
In order to describe the computations required and the accuracies 
obtained, each method is applied to the frequency-response da-ta of fig- 
ure 1 -which defines a single second-order syst.em described by the 
transfer-function 



1 

s^ + 6s + 10 


( 1 ) 


It is usually intended that these methods be applied to hi^er order ' 
systems as -vtIII be shown later in the paper. No attempt is made to 
repeat the development of the methods since this information may be 
obtained from the references. 


Inverse Laplace Transform Methods 

Floyd's method .- The me-thod developed by George F. Flcyd aufl described 
in detail in reference 2 is referred to as Floyd's method. Floyd shows 
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that the Inverse Laplace transform. h(t) of H(s) given lay the 
Integral 


- pc+i» . 

h(t) = — / H(s)e^® ds 

2ltl 


C-loo 


is for all positive values of time equivalent to 


h(t) = / Re |H(ici>)l cos to do 

^ Jn 


QMs operation is based on the assumptions that H(s) may be written 
as the ratio of two rational polynominals in s with real constant . 
coefficients, that lim H(s) = 0, and that H(s) has no poles in the 

S — ^ 00 

ri^t half of the s -plane or on the imaginary axis. The procedure 
for performing the integration required by equation (3) is to plcyfc 
Re[H(io)] against o and then to approximate the exact shape -with a 
series of stral^t-llne segments. The strai^t-line approximation is- 
written as a sum of trapezoidal functions and equation (3) is applied 
to each of the trapezoids; the resulting time functions due to each 
trapezoid are then added to obtain h(t). 

A simple illustration is shown in sketch A: 


![H(ia)^ 



“b 

Sketch A 


“2 % 


The time function associated with a typical trapezoid (as shown by the 
cross-hatched section) may be expressed as 


. * ^ , 
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■where 


A2 = 2^2“^ 


la general for n trapezoids the 


time response is given hy 


li(t) = ^Y~ 



V j 


(5) 


In figure 2 the real part Ee|H(ia3)| of the simple second-order 
system confuted from figure 1 is shown plotted against o>. The func- 
tion h(t) is then easily e'valuated since tables of are given ' 

in reference 2. In the iq)per part of figure 2, two strai^t-line fits 
to the Ee(H(ia)3 aa:e shown, one a five-line fit and the o'ther an eleven- 
line fit. For clarity, however, only the points and not the connecting 
lines of the fit are shown. In the lower part of figure 2 the resulting 
responses to a vmit liipulse coraptcbed by Flpyd's method are compared with 
the exact response to a 'unit Impulse, Ihe accuracy of Floyd ‘s me'thod 
depends on the nimiber of lines used to fit the Ee[s(ia>J] and the loca- 
tion of the cu-toff frequency. !Hie cutoff frequency is defined as the 
maximum frequency at which the Ee[H(ia>2 "was fitted. It should be 
noted that the limits of the definite Integral of equation (5) are 0 
to «•; however. In the practical case the Ee|^H(ia>J] is cut off at 
some finite frequency. 3his error is reflected at the low "values of 
time, especially at t = 0. For a given nuniber of lines, the accuracy 
also depends on the Judicious fit of the lines. - 

ITumerical-integration method . -r la order to use automatic computing 
machines to perform the Inverse Laplace transform method, the necessary 
operations Indicated by equation (5) are performed In the following 
manner by using a numeric8Q.-lntegration method. For a given value of- 
tlme t, the product curve Ee[H(iaj3cos -to is evaluated over the o) 
range nnfl integrated by numerical-integration techniques which give one 
point on the time history of the response to a unit impulse. ly repeating 
the above computation for all the desired values of time, the time response 
can be obtained. 
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This method -was applied to the frequency-response data of figure 1. 
The accuracy of the method is shown in figure 5. In the trErper part of - 
figure 3 the error due to the Interval An chosen for the computation 
is shown by the comparison between the circle and square symbols. The 
accuracy of the computation for the case of Ad = 1.0, as indicated by 
the circle symbols, was not satisfactory beyond t = 2.0 seconds. This 
result is to be expected since the interval Ao> = 1.0 was too large to- 
permit numerical-integration methods to perform adequately the integra- 
tion required by equation (3 ) • 

Rectangilft'r-pulse method .- A method for determining the time response 
to a unit Impulse from frequency-response data has recently been given 
in reference 3- It is referred to herein as the rectangular-pulse method. 
The method involves the iise of tables for time-plane values equivalent; 
to unit rectangul^ pulses of the Be jH(lo>3 . Q3ie method requires fitting 
the Ee[H(lja)j] with a series of rectangular pulses or a staircase func- 
tion so that the area under the curve is equal to the area of the pulse 
in each case as is shown in sketch B: 



Sketch B 


This fitting althou^ made visually should be made carefully. If the 
ordinate of pulse 1 is designated as r^, of pulse 2 as rg, and of 

pulse n as r^, then the time response to a unit Impulse coiresponding 
to this Ee[H(ia))3 is given by 


li(t) = H r.h.(t) (6) 

J=1 ^ ^ 


The time functions h^(t), hg(t), associated with each of 

the rectangular pxilses shown in sketch B are tabulated in reference 3- 
Ohese time functions are the inverse Fourier transforms of vmit rectan- 
gular pulses of the Re|ja:(ia)3 • 

The method was applied to the frequency-response data of figure 1. 
The fitting of the staircase function to the Re[H(iaj)] is illustrated 
in the upper part of figure k and the accuracy of the method is demon- 
strated in the lower part of figure Again the error In the computed 
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re^onae to the unit impulse at t = 0 is due to cutting off the 
Re{H(ia))J at oa = 20 radians/sec and not integrating to <» as indi- 
cated by equation (3) . 


Fourier Method 


Fourier response to imit impulee .- In reference h, the response to 
a \mit step input is derived by the Fourier method. Differentiating 
this response gives an expression by -which the response to a unit impulse 
may be approximated. For a linear system, characterized by its frequency 


response -wl-th amplitude ratio 



and phase angle this expression 


is 




(2n-l)a^ 


cos (2n - l)o^t + 




(7) 


In using equation (7) the choice of the fundamental frequency 

determines the accuracy and length of the computations. Since the 
accuracy is affected by it has been found experience that a 

value ■ 0 ^ = ^ can be \xsed as a ®3od first estimate, -where is the 

undamped natural frequency of the system. Instead of c%, the value of 
frequency at -which the anplittide ratio peaks for a li^tly damped system 
may be used. Ihis -value may be determined from -the frequency response 
of the system. 


In table I the numerical computation of h(t) for the system defined 
by -the frequency response of figure 1 is shown. In this compu-tation, 
ill- terms -were carried in the expression for h(t), and more accuracy may 
be ob-tained, of course, by carrying more terms. Ihe accuracy of -the 
me-thod is shown in figure 5 by -the ccmparison between -the conputed and 
exact response to the unit tnpulse. 


Other Methods 

Schumacher * s method .- Ihe me-thod of reference 5 permits the computa- 
tion of -transfer-function coefficients by assuming -the shape of -the 
■transfer fvmction relating -the input and ou-tput and then curve fitting 
this relation to -the aircraft da-ta in frequency-response form. Once the 
transfer-function coefficients are known, the system is completely 
specified since the response to a unit irpvilse input or any other arbi- 
trary input may then be conpvrbed by -the normal methods a-vailable for 
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solving differential equations. To illustrate the method a san^le 
computation is performed hy using the system defined in figure 1. Ohe 
differential equation relating input 8 and output x for the fre- 
quency response shown in figure 1 is 


+ KjDx + = K58 


( 8 ) 


applying the vector least-squares method of reference 5, the following 
set of simultaneous equations is obtained: 




d=a>L 


- few. 


" feW)^ 




d=u>L 


Oh 








% 




^2 



) (9) 


where 


and 




^ - i Bln 


These equations are then solved simultaneously for the transfer coeffi- 
cients K5, and K^. From, the form of equation (8), it is seen that 

the response to a u ni t Impulse input is given hy the equation 


h(t) 


% 


sin o^t 


(10) 
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■where 


and 



A sample computation Is shown in table U. Table II demonstrates 
the compu-batlonal steps in-volved in sol'ving for the -transfer coeffi- 
cients and for the system, defined in figure 1. These 

transfer coefficients are then substituted into equation (lO), and a 
plot of this function is con^ared with the exact response of the system 
to a unit Inipulse input in figure 6. 

This method ga.ve additional information when applied to the second- 
oider system since the transfer-function coefficients -were determined in 
the process of the confutation. A reasonable amount of ccmputing time 
was- required j and good accuracy -was obtained. The method, however, 
requires pre-vious knowledge of the form of -the transform function 
relating the input and ou'tput. Use of the curve-fitting method on the 
hi^er order system did not prove too effective in determining all the 
parameters of the flexible system. If the form of the transfer function 
of the short period is asapmed and fitted -to ■that portion of the aircraft 
frequency response associated with the short period -then the me-fchod is 
•very effective in determining the short-period transfer-function 
coefficients . 

I*-transform method .- Q3ie P--fcransform me-thod, as described in 
appendix A, differs from the other methods presented in this paper in 
•that the time response to a unit Impulse can be determined directly 
from a known response to a known arbitrary input. The computation is 
carried out entirely in the time domain, bypassing the frequency plane 
entirely, and does not require knowledge of the -transfer function relating 
the input and output. The method also represents a sinple procedvcre for 
using the response .to a unit impu3.se and determining -the response to a 
given arbitrary input. 

An example of the method is shown in figure J in -vdilch the assumed 
input 5(t) and ou-tput x(t) are given from -which the response to the 
unit Impulse h(t) is conptited. A comparison of the computed response 
to a unit impulse and the exact response to a unit Impulse is shown in 
the lower paart of figure J. 



MCA Hff 37Q1 


11 


Response to Arbitrary Inputs 

The procedure for obtaining the response to arbitrajcy inputs -when 
the response to the unit impulse is known is simply a matter of applying 
Duhamel's integral, to the response to the unit impulse and time history 
of the arbitrary Input. A nimerical method of performing this operation 
is shown in appendix A of reference 3. 

The response to a unit step input may be obtained T?y. numerically 
Integrating the response to a unit impulse by using the integrating 
matrix given in reference 3. In like manner, the response to a ranp 
input may be obtained by numerically integrating the rehouse to a step 
function. The' response to a triangular Input is obtained by superposi- 
tion of the responses to various ramp inputs translated along the time 
scale . 

The response to an arbitrary input may be obtained from the response 
to a unit litpulse by the P-transfoim method by use of equation (a 4) of 
the appendix. 

For directly determining the time response to an arbitraj^ Input 
from frequency-response data, the Fourier method is perhaps the best 
known. The general method is indicated In references 6 and 7, In general, 
an Input -vdiich can be eaqjressed as a Fourier series can be represented as 


00 

8(t) “ Aq + ^ ^sln nco^t + 


( 11 ) 


The response to the Input given ly equation (u) of a linear system 
whose amplitude ratio is 1^1 and phase angle is can be written as 




For a few specific shapes of Inputs the values of Bqj Cpci^, and 
of equation (l2) are known. Two wave forms frequently used are the 


square wave and triangular wave and examples of these are given. 
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Response to a square vave .~ As Indicated in reference it the Pourier 
series for a xmit square-wave input 5(t) may he written as 


5(t) = I + a 

2 It 



sln(2n - l)o!pt 


(13) 


and the response to this tmit step input of a linear system of which the 
amplitude ratio is and the phase angle is is given by 


x(t) 





G3ie acctiracy of the confutation d^ends on the choice of the optimum 
fundamental square-wave frequency and, as before, a suitable value is 

usually top = ^ where is the lowest undairped natural frequency of 

5 

the system. 

Response to a triangular-wave Input, — Ihe pacific fom of the input 
used in determining the time response to a triangular-wave input is shown 
in sketch C: 


6(t) 



Sketch C 


where 0^ is the total period of the input, is the length of the 

base in seconds, and 6 is the maximum value of the input. 

Without giving the analytical development, representation of the 
Fourier series for this triangular-wave input is given by 
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8 (t) 


00 


.^0 


2 n-l 



where 


flOp = radlans/sec 

I If 


^8 r , 

3- _ = — V 1 - cos(2ri - l)jt rr- 


and 


^2n-l 



( 2 n - 1 ) 



(15) 


(16) 

(17) 


(18) 


3 he response to this triangular-wave input of a linear system of idiich 
the ampUtude ia ||| ahd the pheae angle la la then 


x(t) 



l)o)^t + 


(19) 


Here again a suitable first approximation is 



For a hi^er order system such as a linear system with several 
structural modes (as Indicated by peaks in the amplitude-ratio curve), 
the choice of a fundamental frequency becomes difficult inasmuch as 

odd mialtiples of cop must give the natural fnequency and the 

frequency of each of the higher structural modes. 
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If axL ipput Is nonperiodic and cannot be expanded into a Fourier 
series directly, tbe Fourier transform of the input is used. Ihe tables 
of reference 5 may be used to determine numerically the Fourier transform 
of an arbitrary input. The proc^iire for computing the response of a 
linear system to such an arbitraiy input is briefly described as follows. 
Ihe Fourier transform of the output is formed by multiplying the Fourier 
transfoonn of the arbitrary Inptit by the frequency response of the system, 
which is the ratio of the Fourier transforms of the output and input. 

Ihe inverse Fourier transform of this product is then taken which is the 
time response of the system to the arbitrary input. Ihe tables of refer- 
ence 3 may also be used to perform the operation of the inverse Fourier 
transform. 


COMPARISON AND USE OF MEJETODS 
Conparlson of Methods 


A comparison of the methods as applied to the second-order system, 
defined by the frequency response shown In figure 1, may be made by 
noting the differences between the exact response to a unit Impulse and 
the response conputed by each of the methods shown in figures 2 to 7* 
itom these figures there ajpears to be little difference between the 
accuracy of any one method over the others, and in each case greater , 
accuracy may be obtained at the cost of more conputation. Fbr the 
comparisons shown, however, the conputlng time required by the 
rectangular-piilse method was significantly less than the time required 
by the other methods. Ihe Fourier response to a miit impulse required 
the most time. A listing of the methods in the order of the conputlng 
time required for the conparisons shown in figures 2 to 7 is given as 
follows; 

Rectangular pulse (ref. 3) 

Schumacher's (ref. 5) 

Floyd's (ref. 2) 

Numerical integration 

Foiirier response to a unit Impulse 

(All computations were performed on a desk-type computer.) ’ 

A more severe test of the methods occurs when they are applied to 
linear systems of higher order. In order to demonstrate an application 
of this lype, a further conparison of the methods was made by applying 
them to the system defined by the transfer function 

x/g\ _ 1 100 

e ' ' p p p 

° s + 6s + 10 s + O.ll-s + 100 s + 0.2s + 225 


( 20 ) 
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The frequency response of this system is shown In figure 8. The real 
part of the frequency response shown in figure 9 was computed from 
figxnre 8 hy using the relation 

Re[H(io>3 = ||cos 03,5 (21) 


In figure 10 a con^arison of the response to a unit impulse computed hy 
four of the methods is shown. For Floyd’s method a 27-line fit to the 
Re[H(ifl))J shown in figure 9 was used’. For the rectangular-pulse method, 
the staircase fit to the Re[H(la>3 'was made with the interval 
Aoj = 1 radian/sec, and for the numerical-integration method, an interval 
An = 0.5 radian/sec was used. She response of the linear system defined 
■fay equation (20) to a triangular input shown in figure 7 of reference 3 was 
used to conq)ute the points shown in figure 10 for the P-transform method. 

A time Interval of At =0.1 second was used. The Fourier response to a 
\mit-inrpulse method was found to he impractical in this case because of the 
large amount of computation required to obtain ajoy accuracy. The accuracies 
of each of the foxn* methods illustrated appear to he equivalent. 


Use of the Response to a TMit impulse as a 
Fli^t-Data-Anaiysis Technique 

As a further comparison of the methods applied to a hi^er order 
system, a typical longitudinal maneuver (M = 0.82) for a flexible sw^t- 
^dng aiiplane has been anaJyzed. In this case the output was the 
pitching-velocity response S at the center of gravity of the airplane 


and the input was the elevator angle 


G3ie analysis of this maneuver 


serves to demonstrate the use of the method of obtaining the time response 
of linear systems to a \mit impiilse from frequency-response data as a 
fli^t-data-analysis technique. The Fourier integrals of the output 0 
and input 6e were evaluated "by using automatic computing equipment and 
the methods of Integration of product curves were used. Hie frequency 
response obtained by dividing the Fourier integral of the output "by the 

Fourier integral of the input is shown in figure 11. !Eie Re|^(icD) 


computed from the frequency response of figure 11 is shown in figure 12. 

In figure 13 a conparison of the response to a unit inpulse conputed by 
three of the methods is shown. For Floyd's method a 31-l.iD.e fit to the 

Rej^(lcD)j shown in figure 12 was used and h^(t) was conputed at enou^ 


values of time to define adequately its shape. For the rectangular- 


pulse method, the staircase fit to the 


was made with the 
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interval ^ = 1 radlan/sec for -the numerical-integration method^ 

an interval of ^ = 0.5 radlan/sec vas used, hi -the case of the 
numerical method^ points are shown In figure 15 only at time Intervals 
of 0.1 second In order to compare accuracies. Several attempts were 
made to compute -bhe response to the unit io^ulse hy -bhe P-transform 
method using -the original elevator input 8e("t) output S(t) daown 

in figure l4. The confutations were made at time intervals ranging from 
0.02 second ^ At ^ 0.1 second hut did not yield satisfactory accuracy. 

It appears that this inaccuracy is primarily due to -fche sensitivity of 
•the me-hhod to aman errors in -the first few terms of -the response or 
Input. 

Time histories of -the original eleva-tor Input Se and pitchlng- 
velocl-ty response 0 are shown in figure llj-. The eleva-tor input Be 
has also heen vised as a forcing function wl-th -the response -to a unit 
Ijifulse shown in figure 15 -to conpute a time history of 9. The Duhamel 
me-thod outlined in appendix A of reference 5 wi-th a time Interval of 
0.05 second was used for -these calculations. A conparlson of -this 
confuted 0 response -wi-th -the original 0 response from a fli^t 
test is shown in -the lower part of figure 14- . This comparison shows -the 
amount of error Involved in -the to-tal computation procedure (-transferring 
the da-ta In tlme-his-tory form to frequency-response form, -then to -the 
response to a unit Inpulse, and finally to -the response to an arhi-trary 
input hy means of Duhamel* s integral), 

Ano-ther fli^t-test maneuver at M = 0.80 -was selected -wi-th -the 
o-ther conditions approximately -the same as -the previous da-ta In order 
to see how well -the response -to a unit lufulse computed from one maneuver 
could he used -to predict -the time response -to an arhi-trary input from a 
different maneuver. The time his-tories of this maneuver are shown in fig- 
ure 15. Also shown in -the lower part of figure I5 is -the 5(t) response 
calculated hy application of Dvihamel's integral to -the response to a unit 
Impulse given in figure 15, and -the elevator motion for the maneuver is 
shown In -the xfper part of figure 15. A confarison of the computed and 
measured fi response in -the lower part of figure I5 gives some indication 
of how well -the motions of an aircraft can he predicted hy a de-talled 
analysis of a single maneuver. 


COBCLUPING RE61AEKS 


Se-veral methods -were compared for oh-tainlng -the time response of 
linear systems -to ei-ther a unit la^julse or arhi-trary input from frequency- 
response da-ta. The methods -were compared on -the basis of accuraoy, 
conputing time required, and applicahill-ty to higher order linear systems. 
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The application of each of the methods to a simple second-order 
system indicated little difference between the accuracy of one method 
over the others, and, in general, it would be expected that greater 
acciaracy could be obtained for each of the methods at the cost of more 
confuting time . 

For higher order systems the three methods genereilly classified as 
inverse Laplace transform methods were most effective. They gave good 
results for a moderate amount of computation. All of these methods are 
based on the evaluation of the form of the inverse Laplace transform 
equation for positive values of time. 

The methods generally classified as Fourier methods gave good 
accuracy when applied to the second-oirder system. The accuracy of these 
methods was found to depend on the choice of a ftmdamental frequency o^. 

For simple systems, a suitable value was found to be ^ where 

is the lowest natural frequency. For higher order systems the choice 
of oif becomes more difficult. The Fom:*ier response to a unit-impulse 
method was found to be in^iraoticaJ. in the case of hi^er order systems 
because of the large amount of computation required to obtain accuracy. 

Schumacher's method gave additional Information -vdien applied to the 
second-order system since the transfer-function coefficients were deter- 
mined in the process of the computation. This method required a reason- 
able amount of computing time and gave good accuratqr. The method, 
however, requires previovis knowledge of the fom of the transfer func- 
tion relating the input and output. Use of the method on higher order 
systems did not prove too effective in determlng eJd the parameters of 
the flexible system. A valuable use of the method was found in fitting 
the known short-period transfer function to only the short-period portion 
of the flexible-system frequency response to determine the short-period 
transfer coefficients. 

The P-transfoim method is different fi^xii the other methods presented 
in that the time response to a unit impulse can be determined directly 
from a known response to a inown arbitrary input. The computation is 
carried out entirely in the time domain and bypasses the frequency plane 
entirely. The method also r^resents a simple procedure for using the 
response to a unit Impulse and determining the response to a given 
arbitrary input. The method, however, when applied to the fli^t data 
of this paper did not yield satisfactory accuracy. It appears this is 
primarily due to the sensitivity of the method to errors in the first 
few terms of the response or input. 

For all the examples conputed in this paper the rectangular-pulse 
method of MCA Technical Uote 3598 required less computing time, for the 
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same degree of accuracy, -than aaay of the other methods used to obtala 
the time response to a unit ln^julse from frequency-response data. 

Men an aircraft frequency response Is known, it aj^ears that these 
methods may be used as a fligit-data-analysis technique which permits 
prediction of aircraft motions and loads without knowledge of the equa- 
tions of motion relating the inputs and outputs for a flexible aircraft. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Tangifiy Field, Va., March 16, 1956 . 
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APPENDIX A 


P-TRANSK)BM METHOD 


In reference 8. a linear operational calculus is Introduced -which 
appears to be well aidapted to the numerical analysis and synthesis of 
linear systems. In this calculus a polynomlnal transform or P-transform 
of a f motion F(t) is defined by the equation 

PlF(t)l =f(x) 

= P(mAt)x?^ (Al) 

m=-K 

where P(mt) are the ordinates of P(t) at integral multiples of a 
time InterWl At. The inverse P-transform is given by 

F(t) = P"1 [f(x3 (A2) 


The si^erposition (convolution co* Ehltung) integral, which is also known 
as Duhamel's Integral, relates -the Input 6(t), output x(t), and 
response to a mit trapiilse h(t) (also called a memory ftmc-fcion) of a 
linear system and is given by 


x(t) = f 8(t) h(t-T)dT (A3) 


Bubb (ref. 8) shows that the P-transform of equation (A3) Is 

x(x) = At b(x) h(x) (Alf) 


In the synthesis problaa the Input 6(t) and the ou-i 5 >ut x(t) are 
given and the response to the unit impulse h(t) is to be ceilculated. 
This calCTilatlon is perfoimed by forming the P-transfoims of x(:^ and 
&(x) and dividing by ordinary polynomlnal di-vLsion, x(x) by 8(x), to 
get the P-transform of the response to the unit Impulse 
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h(x) 


For practical cases, however, _ 

of t values. Tfcie inverse P-transform to h.(x) then gives h(t), the 
time response to the unit imptilse. An example of this operation is shown 
in table TTT ; the known input 8(t) and output x(t) are shown in the 
upper part of figure 7^ a°d the confuted response to the unit impulse is 
comgpared with the analytic solution in the lower part of figure 7* 

Thifi operational calculus is also well adapted to the solution of 
the analysis problem in which the system response to a unit impulse h(t) 
and the input function S(t) are known and the calculation of the system 
output x(t) is desired. The P-transforms of h(t) and 5(t) are 
formed anfl multiplied together by ordinary polynominal multiplication as 
Indicated in equation (A^-) . Since this operation is Just the Inverse of 
the operation shown in table TTT , an illustrative example is not shown. 
This method has been found to be a simple and rapid means of applying 
the Duhamel process. 

The value of this operational calculus lies in the fact that all 
computations remain In the time domain and no translation to the frequency 
plane is required. Also, only sln^jle direct arithmetical procedures are 
required for solving practical problems. 


1 x(x) 
At 5(x) 




0 “ 


(A5) 


he summation is made over a f Inite range 
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TABLE I 

HIMERICAL CALCULATION ILLUSTRATING COMPUTATION OP THE 
FUNCTION FOR THE RESP05BE TO A UNIT IMPULSE 
BT THE FOURIER METHOD 


0^ = j 0^ = 0.6^ 


n 

2n - 1 

(2q - l)a^ 

lx 


^8^ 

deg 

1 

1 

0.64 

0.0970 

-22 

2 

3 

1.92 

.0760 

-61 

3 

5 

3.20 

.0520 

-91 

4 

7 

4.48 

.0350 

-m 

5 

9 

5.76 

.0240 

-124 

6 

11 

7.04 

.0170 

-153 

7 

13 

8.32 

.0130 

-l40 

8 

15 

9.60 

.0100 

-145 

9 

17 

10.89 

.0079 

-149 

10 

19 

12.16 

.0064 

-152 

11 

21 

13.45 

.0053 

-155 

22 

23 

14.72 

.0044 

-157 

13 

25 

16.00 

.0038 

-159 

14 

27 

17.30 

.0033 

-160 


h(t) = ^^^[0.097 cos(0.61tt - 22) + 0.076 cos(l.92t - 61) + 
oc 

0.052 cos(3.2t - 91) + 0.035 cos(li-.48t - 111) + 

0.02l^ cos(5.76t - 124) + 0.017 cos(7.04t - 135) + 
0.013 cos(8.32t - l40) + 0.010 cos(9.6t - 145) + 
0.0079 cos(l0.89t - l49) + 0.0064 cos(i2.l6t - 152) + 
0.0053 cos(l3.45t - 155) + 0,0044 cos(l4.72t - 157) + 
0.0038 cos(l6t - 159) + 0.0033 cos(l7.3t -160)] 
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0.1000 

.09^ 

• 07^^5 

•0555 

.o4oii- 

.0298 

.0225 

.0174 

.0139 



-104.0 

-116.6 

-125.8 

-332.9 

-138.4 


0.10000 

.07687 

.03336 

.00310 

-.00977 

-.03334 

-.01316 

-.01184 

-.01039 


^ “ II sin 


0 

-.05127 

-.06661 

-.05541 

-.(0920 

-, 6 s 6 &v 

-.01825 

-.01274 

-.00923 


- 






^ /|x|2\ 


K5 = 1.0018 
Ki = 6.01865 
= 10.0072 


-0.913^ 

-0.15483 


-0.933^ -0.15483 % 


O.152OL 


O;C0O69 


1.97049 


0.15201 


li(t) = ^ e"®* sin o^t = 1.03e“^’°°^^ sin 0.9754 


= jKe - 


a 






mu m 


mHERBUL oALouuzioir iLLDssBAxan me t^RAsaFoai uEmiD 


[*'“> ‘il§j 





Oul^^, 

3t(x) 


.oooaJte. 

.001X0*8 

.00292*3 

.ooeos*?*- 

.01051x5 

.0153TK® 

.01994:J 

,02516*0 

.02^66*9 

.0*33x^? 

,cKago*“ 

.01681x?5 

.015933^ 

.0U3O»W 

.009®*^ 



/ . . o ^ at k . — ^ 

^ 1 J O + o«qootfx + o.ooioor* » o^ooggac^ + g«womjx^ -h u*{uuw' * • • » + u« 
0.l\ 0 ■!• O.lx + O.Ei^ + 0.3a? + 0.4s^ + . . . + 0.1x9 

E(x) » O.0o;i' + O.OTSx -I- O.10&^ *■ 0 . 1191 ? + O.UT*^ ■•■...- 0.090a^ 


















Frequency ,:w , radians/sec 


Figure 1.- Frequency response of second-order system, 


1 

s^ + 6s + 10 





Response to a unit impulse, h(t) Re[H(icjj| 


26 


MCA TN 5701 



Figure 2 .- Time response to a vmlt impulse obtained from Floyd's method 
compared to the exact values for the second-order system. 




Response too unit impulse,h(t) Re[H(icj8 
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jy_gure 4 .- Hme response ‘to a vinil; Impulse obtained from the rect an g ul ar - 
pulse method compared to the exact values for the second-order system. 


% 7 9 





Figure Time response to a unit lit 3 >^LLse obtained from the Fourier 
method conqpared to the exact values for the second-order system. 
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Figure 8.- Frequency response of higher order system. 

X/g) ^ 1 100 225 

S s^ + 6s + 10 s^ + 0.4s + 100 s^ + 0.2s + 225 


Re[H(i(u)] 
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Figure 10 .- SPime response to a unit impulse obtained lay four methods 
CQo^ared to the eacsuit values for the higher order system. 
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Figtire 11.- A typical frequency response obtained from fli^t data of a 

flexible airplane at M = 0.82. 


Re [4 (icu}l , radians/sec/deg 
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Figure 12.- Real part of the typical frequency response obtained from 
flight data of a flexible airplane at M ■> 0.82. 
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Pitching-velocity response to a unit impulse. 
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Figure IJ.- Time reaponae to a unit impulse obtained Isy three methods 
for the typical frequency response obtained from flight data of a 
flexible airplane at M 0.82. 


I 


mCA TN 5701 



58 


NACA TN 3701 




Figure l 4 .- Coraparison witli the actual flight record of the calculated 
pitching-velocity response using the same iiq>ut-ou-^ut data as was 
used to obtain the typical frequency response of the flexible air- 
plane at M = 0.82. 
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Figure I5.- Camparlson of the predicted pitching-velocity response 
with the actual flight record obtained from the flexible air- 
plane at M = 0.80. 
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